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Abstract. Both, X-ray and radio observations, are beginning to probe regions in galaxy 
clusters close to the virial radius. These are particularly interesting regions as they are (i) 
crucial for the thermal evolution of the intra-cluster medium and (ii) they host plasmas 
under conditions that are poorly understood. Recently discovered giant radio relics in the 
galaxy cluster CIZA J2242.8+5301 provide unprecedented evidence for particle accelera- 
tion in cluster outskirts. We will discuss some of the intriguing puzzles that surround these 
relics. The polarisation of the radio emission yields information about the topology of mag- 
netic fields in these shock waves. Here we report on MHD simulations that highlight some 
problems about the origin of these magnetic fields. 
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1. Introduction 

Radio relics are diffuse radio sources that are 
usually located in cluster outskirts and that 
are characterized by steep radio spectraQ with 
a > 1.2. The origin of the relativistic electrons 
- that produce the synchrotron radiation in the 
radio band - in these very extended sources is 
still not entirely clear. In the past two years, 
there has been quite a bit of new data on ra- 
dio relics with the discovery of the relic in the 
galaxy cluster CIZA J2242.8+5301 being the 
most prominent. CIZA J2242.8+5301 is a dis- 
turbed galaxy cluster located at z - 0.1921 
dKocevski et al.l 120071) . The cluster hosts a 
large double radio relic system, as well as 
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additional large-scale diffuse radio emission 
(Ivan Weeren et al.l [20101) . This recent data has 
given support to models that explain relics 
through diffus ive particle acceleration at shock 
fronts , e.g. (lEnsslin et al.l Il998b iMiniati et all 
2001). The spectral index (a) across the bright 
northern relic steepens systemically in the di- 
rection of the cluster center, across the full 
length of the narrow relic. This is expected 
for outwards moving shock waves, with syn- 
chrotron and inverse Compton losses behind 
the shock front. However, there are still quite a 
few puzzles that surround radio relics, of which 
I will list the most intriguing: 



- Which processes accelerate electrons so ef- 
ficiently at relatively low Mach number 
(M ~ 2 - 4) shocks? 
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- What produces the magnetic fields inside 
relics? Both the inferred field strengths as 
well as the observed polarisation of the ra- 
dio emission need to be explained. 

- Why do some relics have very sharp edges 
while others appear very fuzzy? 

- Under which conditions do relics form? 
When do we see single and when double 
relics? 

- Some relics appear to be connected to 
cluster-wide radio halo emission. Is there 
a physical connection between the two? 

Many of these questions touch physical 
processes that are poorly understood, such as 
diffusive shock acceleration and magnetogen- 
esis in collisionless shocks, and relate to re- 
gions of the cosmos, namely clusters outskirts, 
of which we know very little. Ongoing work 
with radio interferometers, in particular with 
the new LOFAR radio telescope is likely to 
bring great advances in these areas. 

In these proceedings, we are going to 
present some recent theoretical work that ad- 
dresses the orientation of magnetic fields in 
cluster outskirts. 



2. Simulations of mergers 

Idealised simulations of cluster s hav e 
been presented by FRoettieer et al.l (1 1 9931) : 
ISchindler & MueUerl dl993b : IPearce et al.l 



(1994]); IRoettiger et~ai1 dl997l) . iRoettiger et al.l 
(1 19931) simulated the head-on collision of 
galaxy clusters with unequal mass, each of 
them modelled as King- spheres. They used 
the ZEUS code with a non-uniform grid and 
the Hernquist (1987) treecode for the dark 
matter. The simulations showed a bar of 
X-ray emission perpendicular to the collision 
axis that coincided with a shock of M ~ 4. 
Schindler & Muellei] (Il993l) used a uniform 
grid-PPM code + an N-body code to study 
the head-on collision of galaxy clusters. 
They observed several shocks that lead to a 
substantial enhancement in X-ray luminosity 
(up to a factor of 2). Their simu lation did 
not in clude cooling. Later on, iBurns et aD 
(1994) simulated the collision between the 
Coma cluster and the group NGC 4839, using 



the ZEUS code with radiative cooling, and 
obtained similar results. 

None of these works included the baryonic 
contribution to the gravitational potential. 

IPearce et aD (1 19941) simulated the head- 
on collision between two equal-mass clus- 
ters using smoothed-particle hydrodynamics 
(SPH) that includes, both, baryons and dark 
matter. They found that the baryons end up 
in extended constant density cores after the 
collision, wher e as the dark matter does not. 
IRoettiger et~ai1 (1 19971) used a PPM/particle- 
mesh code to simulate the radio relic in A754, 
varying the impact parame t er and mass ra- 
tios of the collision. iRicked (1 19981) presented 
simulations of off-centre collisions between 
clusters using a PPM code with an isolated 
multigrid solver for the gravitational poten- 
tial. They studied how the virialisation time, 
the structure of the merger and the X-ray 
luminosity enhancement depends on the im- 
pact parameter of the collision. They find 
that the X-ray luminosity increases by a large 
factor, which strongly depends on the im- 
pact parameter, and lasts for about half a 
sound crossing time after the cores inter- 
act. Virial isation takes abou t 5 so und cross- 
ing times. iRicker & Sarazinl (l200ll) presented 
a parameter study of off-centre collisions be- 
tween galaxy clusters of two different mass 
ratios. Using the Eulerian hydrodynamics/N- 
body COSMOS code, they collide hydrostatic 
clusters with a ft model gas profile. They note 
that the jumps in X-ray luminosity and tem- 
perature have important consequences for cos- 
mological parameters derived from the ob- 
servations of ho t clu sters at high redshift. 
Rick er & Sarazinl (l200ll) also find that merg- 
ers lead to large-scale turbulent motions with 
eddy sizes of several hundred kpc and this tur- 
bulence is driven by dark matter-driven oscil- 
lations in the gravitational potential. 



2.1. Method 

The simulations were performed with the mul- 
tidimensional adaptive mesh refine ment hy- 
drody namics code FLASH code (iDubey et al.l 
2009) version 3.2. 
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A directionally unsplit staggered mesh al- 
gorithm (USM) that solves ideal and non-ideal 
MHD governing equations in multiple dimen- 
sions is newly implemented in this latest ver- 
sion of the code. The overall procedure of the 
USM scheme can be broken up is described in 
the FLASH manual and in Lee (2006) and Lee 
and Deane (2007). 

We model the ICM as ideal gas in hydro- 
static equilibrium. The density follows a ft pro- 
file 

p(r)=p [l + (r/r ) 2 r W2 , (D 

where ft = 0.3, po = 10" 25 g/cm 3 and ro = 70 
kpc are parameters. 

The temperature profile is taken to be 

T(r) = (r +ar)[l-(l+exp(-(r-r,)/^))- 1 ],(2) 

where the parameters have been chosen to be 
T = 3 x 10 7 K, a = -5.6 K/pc, b = 10 kpc and 
r\) - 27 kpc. 

The gravitational potential follows from 
the assumption of hydrostatic equilibrium. 

dr jutih \ p dr dr) 

We move the centre of the gravitational po- 
tential of cluster 1 on a trajectory around the 
fixed potential of cluster 2 and thus ignore all 
interactions between dark matter. 

Furthermore, the hydrogen number density 
was assumed to be related to the electron num- 
ber density as nn = n Q /1.2. We did not include 
the effects of radiative cooling or star forma- 
tion, as these are not very relevant for the spa- 
tial and temporal time scales considered here. 

We set up an initially tangled magnetic 
field that follows the thermal component ac- 
cording to: 

(B)(r) = (Bo>(^)" (4) 

where (Bo) is the mean magnetic field strength 
at the cluster center and no is the number den- 
sity at the cluster centre. 

Simulations start considering a power 
spectrum for the vector potential A in the 
Fourier domain: 

\A k \ 2 oc k~{ (5) 



and extract random values of its amplitude A 
and phase 0. A is randomly extracted from 
a Rayleigh distribution (in order to obtain a 
Gaussian distribution for the real magnetic 
field components, as frequently observed), 
while varies randomly from to 2n. All three 
components of magnetic field were treated in- 
dependently which ensures that the final distri- 
bution of B(r) has random phase. 

The magnetic field components in the 
Fourier space are then obtained by: 

B(k) = ikxA(k). (6) 

The field components B ( in the real space 
are then derived using a 3D Fast Fourier 
Transform. The resulting magnetic field has the 
following properties: 
1) 

V • B = 0, 

2) the magnetic field energy density associated 
with each component B k is: 

\B k \ 2 = C\k-\ 

n = £ - 2, where C 2 n is the power spectrum 
normalization, 

3) B t has a Gaussian distribution, with {B t ) = 0, 
cr Bi = {B% 

4) B has a Maxwellian distribution. 

For n = 1 1 / 3, the above prescription would 
give a Kolmogorov distribution of magnetic 
fields. 

One can define A = y as the physical scale 
of the magnetic field fluctuations in the real 
space. Thus in order to determine the magnetic 
field power spectrum in the cluster, we have to 
determine three parameters: A min , A max and n. 

There are several indications that the 
magnetic field intensity decreases going from 
the center to the periphery of a cluster. This 
is expected by magneto-hydrodynamical 
simulations (see e.g. Dolag et al. 2008) and 
by spatial correlations found in some clusters 
between thermal and non-thermal energy 
densities (Govoni et al. 2001). 
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Fig. 1. Initial conditions: slices of logarithmic 
gas densities, magnetic fields, plasma J3 and 
pressure. 

3. Results 

Here we only have space to show a small part 
of our results. For further results we refer to 
an upcoming publication (Bruggen, Lee, van 
Weeren, Rottgering in prep.). In Fig. 1 we show 
our initial conditions and in Fig. 2 we show 
slices of the gas density and magnetic field vec- 
tors in a head-on collision of two galaxy clus- 
ters at various times. These are quite represen- 
tative of many of our simulations. Our main 
findings can be summarised as follows: 

- Initially the magnetic fields become fila- 
mentary in wake of galaxy clusters as result 
of stretching. 

- As the clusters approach each other: en- 
hancement of density, temperature and 
magnetic field in the midplane magnetic 
field is perpendicular to the shock axis 

- After centre-passage: shock waves run dia- 
metrically out of the cluster. Magnetic field 
is perpendicular to shock front. 

- Cluster gas dynamics gets dominated by 
large-scale flows and large eddies that 
shape the morphology of the cluster-wide 
field. Along the merger axis, the flow is 
largely directed outwards (away from the 
gravitational centres) along the axis. This 



Fig. 2. Slices of logarithmic gas densities (in 
a zoomed region) with magnetic field vectors 
superimposed for an equal-mass head-on colli- 
sion. 

results in magnetic fields that are per- 
pendicular to the outgoing merger shock 
waves. 

- As the magnetic fields observed in the giant 
radio relic in CIZA J2242.8+530 lie within 
the shock plane, these magnetic fields are 
difficult to explain by a mere amplification 
of ambient magnetic fields. 

These findings are consistent with what 
has been found in MHD simulations by 
dRuszkowski et al.l 120101) . They have per- 
formed full cosmological simulations of 
galaxy cluster formation that simultaneously 
include magnetic fields, radiative cooling and 
anisotropic thermal conduction. The motiva- 
tion behind this work was that in isolated and 
idealized cluster models, the magnetothermal 
instability (MTI) tends to reorient the mag- 
netic fields radially whenever the tempera- 
ture gradient points in the direction opposite 
to gravitational acceleration. Using cosmolog- 
ical simulations of the Santa Barbara cluster 
Ruszkowski. et al. (2010) have detected a ra- 
dial bias in the velocity and magnetic fields. 
However, it was found that such radial bias 
is caused by the accretion of gas and sub- 
structure and not by MTI-driven field rear- 
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rangements. The anisotropy effect is poten- 
tially detectable via radio polarization mea- 
surements with LOFAR and Square Kilometer 
Array and future X-ray spectroscopic stud- 
ies with the International X-ray Observatory. 
Furthermore it was demonstrated that radiative 
cooling boosts the amplification of the mag- 
netic field by about two orders of magnitude 
beyond what is expected in the non-radiative 
cases. This effect is caused by the compres- 
sion of the gas and frozen-in magnetic field 
as it accumulates in the cluster center. At z = 
the field is amplified by a factor of about 
10 6 compared to the uniform magnetic field 
evolved due to the universal expansion alone. 
Interestingly, the runs that include both radia- 
tive cooling and anisotropic thermal conduc- 
tion exhibit stronger magnetic field amplifica- 
tion than purely radiative runs, especially at 
the off-center locations. In these runs, shal- 
low temperature gradients away from the clus- 
ter center make the ICM neutrally buoyant. 
Thus, the ICM is more easily mixed in these re- 
gions and the winding up of the frozen-in mag- 
netic field is more efficient resulting in stronger 
magnetic field amplification. 
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